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Abstract 

The analytical description of the volume reflection of the charged ultrarelativistic 
particles in bent single crystals is considered. The relation describing the angle of 
volume reflection as a function of the transversal energy is obtained. The different 
angle distributions of the scattered protons in the single crystals are found. Results of 
calculations for 400 GeV protons scattered by the silicon single crystal are presented. 

1 Introduction 

The volume reflection of the charged particles in single crystals represents the coherent scat- 
tering of these particles by planar or axial electric fields of bent crystallographic structures. 
For the first time, this effect was predicted in Refs.[l, 2]. The angle of the volume reflection is 
small enough (in order of the critical angle of the channeling). Because of this, it is thought 
that this process did not show itself in experimental measurements. However, in the recent 
reports [3, 4] the problems were discussed, which require the existing of volume reflection 
effect for explanation of proton collimation measurements. Besides, the preliminary results 
of the experimental observation of the volume reflection were published [5, 6]. 

The recent meeting in CERN [7] which was devoted to problems of utilization of interac- 
tions of proton beams with single crystals attract considerable interest to the phenomenon of 
volume reflection due to new possibilities such as the collimation of collider beams [8, 9, 10], 
creation of devises which can multiple the effect [11] and so on. In the immediate future the 
experimental study of the volume reflection effect are scheduled on extracted proton beams 
in CERN [12] and IHEP[13]. 

By now the main results of the theoretical investigation of the volume reflection are mainly 
based on the Monte Carlo calculations. No doubt, this is very powerful method, however, 
the analytical solution of the problem is very useful for understanding of the process and its 
different characteristics.. In this paper we consider such description of the volume reflection 
effect . 

2 Motion in bent single crystals 

The motion of the ultrarelativistic particles in the bent single crystals one can describe with 
the help of the following equations[14, 15] : 

E /3 2 v 2 /(2c 2 ) + U(r) + E /3 2 (r - R)/r = E = const (1) 

dy I dt = v y = const (2) 
v z = rd^/dt « c(l - (l/ 7 2 + v 2 + v 2 )/2) (3) 
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These equations take place for cylindrical coordinate system (r,(f),y). Here v r is the com- 
ponent of particle velocity in the radial direction, v y is the component of the velocity along 
y-axis and v z is the tangential component of the velocity, R is the radius of bending of the 
single crystal, E and 7 are the particle energy and its Lorentz factor, E is the constant 
value of the radial energy, U (r) is the one dimensional potential of the single crystal, c is the 
velocity of light and (5 is the ratio of the particle velocity to velocity of light. In this paper 
we consider the planar case when the scattering is due to interaction particles with the set 
of the crystallographic planes located normally to (r, 0)-plane . On the practice it means 
that Vy/c ^> 9 ac but v y /c ^> 1 for ultrarelativistic particles, where 9 ac is the critical angle of 
the axial channeling. 

Eq.(l) one can transform in the following form: 

E (3 2 v 2 J{2c 2 ) + U(x) + E [3 2 x/R = E (4) 

where x is the local Cartesian coordinate which connected with the cylindrical coordinate r 
through the relation x = r — R and v x — v r . We also change r- value in the denominator of 
Eq.(l) on R. For real experimental situation it brings negligible mistake (in order of x/R). 
In Eq.(4) E- value have a sense of the transversal energy. 

The considered here equations describe the three dimensional motion of the particle in 
the bent single crystal in the cylindrical coordinate system. Let us introduce the Cartesian 
coordinate system in which the xy-plane is coincident with the front edge of a single crystal. 
Now we can calculate x-component of the velocity in this system: 

Vx = v x cos + v z sin ss v x + v z (f> (5) 

Fig.l illustrates the geometry of proton beam volume reflection and different coordinate 
systems which will be used in our consideration. 



3 Angle of volume reflection 

One can see that Eq.(4) describes one dimensional motion in the effective potential U e (x) = 
U (x)+E /3 2 x/R. There are two different kinds of the infinite motion in the effective potential. 
One case corresponds to motion when particle after its entrance in single crystal move in 
the direction of the increasing of potential ( in positive direction of x-axis). In this case the 
velocity v x decrease and in some point x c become equal to zero and then particle begin to 
move in the opposite direction. It corresponds to case when effect of volume reflection is 
possible. In the second case the particle after entrance in single crystal move in the negative 
direction of the x-axis and volume reflection is impossible. Obviously, that the belonging 
of a trajectory to some kind of events is determined by the initial conditions of particles at 
entrance in the single crystal: the first case takes place at the negative entering angles 9 
relative to crystallographic planes. We also will consider below case when the absolute value 
of 9 is large than channeling angle 9 C . 

The trajectory in x-direction (xo < x < x c ) up to point x c one can calculate as 

r x dx 
t= / ; o 6 
J*o JM,(E - U(x) - E n 3 2 x/R) 



L ^(E-U{x)-E Px/R) 

where t is the time and t — corresponds to coordinate x at entrance of the single crystal. 
We can rewrite this equation in the following equivalent form: 

* = 1 = - F(x)]dx + f F(x)dx (7) 

Jxo [ ^ {E -U(x)-E ^x/R) J*o 
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where the function F(x) reads 



F(x) = . , = (8) 

^(E-U(x c )-E ^x/R) 

and the critical point x c is determined by the equation: 

E - U(x c ) - E (3 2 x c /R = (9) 

Note that equation t = f* F(x)dx describes the one dimension motion of a particle, which 
do not interact with the electric field of the single crystal. We name this case as the neutral 
motion. 

Taking the integral f* F(x)dx we get 

t = T(x ,x) + R/c 2 (v(x )-v(x)) (10) 

where 

T(x ,x)= /*[ i o 1 =-F{x)}dx (11) 

1 ' 2c2 'E - U{x) - E f3*x/R) 



and v(x) is the velocity of the neutral motion in the point x of its trajectory: 

v(x) = ^2c 2 /(E ^)(E - U(x c ) - E (Px/R) (12) 

The angle in Eq.(5) is approximately equal to V z t/R ps Vt/R ps ct/R For our aims the 
small difference between V z ,V,c is of no significance. We will use the relation: = ct/R. 
Hence, Eq.(5) have the following form: 

V x = v x + v z ct/R&v x + c 2 t/R (13) 

Taking t from Eq.(10) one can get 

V x = v x + T(x , x)c 2 /R + v(x ) - v(x) (14) 

Now we assume that the distance d vr — \x c — Xq J is large enough. According to estimation 
in Ref.[2] d vr ^> 6 C R. Hence, V x ~ v(x ) (v x ~ v(x), T(xq,x) ~ 0) when the particle move 
far from the critical point (near x$ point). It easy to see that in the critical point v x = 0, 
v(x) = and V x = v(x ) + T(x ,x c )c 2 /R. Using the obvious symmetrical properties of the 
one dimensional motion we can get that for distant enough from critical point values x > x c 
V x = v(x ) + 2T(x , x c )c 2 /R. From here the angle of volume reflection a = 2T(x , x c )c/ R or 
its explicit value have the following form: 



a = 4: Hi ; , - F(x)]dx (15) 

RJxo \rM7(E -TT(<A - E„Wt/R) V n V ' 



2c 

' xo ^^(E-U(x)-E ^x/R) 

It is easy to understand the meaning of the F-function introduced in Eq.(8) . This 
function allows one to eliminate the influence of the geometrical factors on the particle 
motion in the cylindrical coordinate system. The visual picture of volume reflection may be 
also observed in cylindrical coordinates r and 0. Really the trajectories of the charged and 
neutral particles have the same critical radial coordinate r c but the different values of C . 
The difference of c -coordinates determines the angle of volume reflection. 
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4 Periodicity of volume reflection angle 

Considering the volume reflection angle as a function of the transversal energy E one can 
see that this function is periodic. The period is equal to 

5E = E (3 2 d/R. (16) 

This fact one can get by the substitution x = x + d (d is the interplanar distance) and taking 
into account the periodicity of potential U (x) . 

The transversal energy E have the following representation: 

E = E (3 2 9 2 /2 + U(x ) + E (3 2 x /R (17) 

The replacement x on x + d corresponds to new energy E + 5E. It means that the angle 
a is periodic function of the initial coordinate xq. 

The volume reflection angle is also a periodic function of the square of initial angle 9. In 
this case the period is 

59 2 = 2d/R (18) 

It should be noted that periodicity of volume reflection angle relative to above pointed param- 
eters violated if the condition d vr ^> 9 C R is not satisfied or, in other words, at transformation 
from one energy to another the distance \x — x c \ should be large enough. 



5 Particle distributions at volume reflection 

In calculations we will assume that the distribution of the particle beam over the initial 
coordinate Xq is uniform (in the limits of one period). If the function a(xo) (0 < x < d) is 
known we can find the distribution of particles over the angle of volume reflection with the 
help of the following relation: 

da y ' ^ dUxo y ' 

where the sum over i means that the sum should be taken over every domain of the function 
single- valuedness and the derivative da/dxo should be calculated for values XQi which satisfy 
to equation a(xoi) = a in which a is the current (computed) value of the volume reflection 
angle. Note Eq.(19) was obtained at the condition that the initial angle 9 is fixed. Besides, 
in this paper we use the normalized to 1 the distribution functions. 

In actual practice, the particle beam have some distribution over the initial angles. Ac- 
cording to our consideration the distribution function over a dependent on the initial angle 
value. On the other hand, the volume reflection angle is the periodic function of # 2 -value. At 
large enough ^-angles and small range of variation of this angle the angle a one can consider 
as periodic function of the initial angle 9 with the period: 

se =m < 20 > 

It easy to see that 89 is very small value. Simple estimations give values of 59 ~ 10~ 6 — 10~ 8 
radian for radius in the range 1- 100 m, correspondingly, for initial angle about 10 -4 radian 
and d about some angstroms. The real particle beams have significantly more a broad angle 
distribution. Assuming the uniform distribution of the initial angles within the angle period 
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89 (at the assumption that divergence of beam satisfied the condition A9 <C 6) we get the 
distribution averaged over its period 

1 r 59 dN 

< dN/da >= — / =—d0 (21) 
89 Jo da 

This equation is obtained for physically narrow particle beam. Under this term we under- 
stand the beam with the small divergence of initial angles, which in tens or more times 
exceed the angle period in Eq.(20). The distribution of the beam with a valuable quantity 
A9 (however, A9 -C 9) have the following form 

dN /-02-0O dN 

where Pb{£,,9o) is the angle distribution of beam around the initial angle 9q, 8i and 9 2 are 
minimal and maximal values of the initial angles (A9 = 9 2 — In this equation the angle 
a is the planar angle relative to direction which is determined by 9 . 



6 Influence of multiple scattering 

The above-mentioned results of our consideration are appropriate for the ideal case only 
when the dissipative processes in the body of single crystal are absent. The Monte Carlo 
calculations can give detail consideration of these processes. Here we make attempt to take 
into account the influence of the multiple scattering on the process of volume reflection in 
the frame of the analytical description. However, in this case some realistic assumptions is 
needed. As was said earlier the main contribution in integral describing the angle motion 
of particle (see Eq. (15) ) give relatively not large a space area near the critical point x c . 
Hence, we can assume that we can neglect by the myltyple scattering in this area. 

The process of multiple scattering on the particle path from xq till ~ x c bring the variation 
of the transversal energy E. Bearing in mind that the scattering can play the important role 
in the small area near x c -coordinate we suppose that this factor may be taking into account 
by the corresponding changing initial data (see Eq.(17)). We assume that the physically 
narrow beam have the initial angle 9. Then we can consider that the initial distribution of 
particles is following: 

dN f°° dN 

<—(a)> ms = p(ip,a) < —(a - ip) > dip (23) 
da J~oo da 

where p{ip) = ^7=^ exp — (p 2 /2a 2 ) is angle distribution of multiple scattering with the pa- 
rameter o depending on the distance ~ \x c — xq\/9. We also expect that the redistribution 
of the initial coordinates is small and do not consider this factor. Here it is important that 
the distribution of the particles within interplanar distance remains uniform. 

Now we can obtain the angle distribution after single crystal for physically narrow beam: 

dN f°° dN 

< —{a) > exit = / p(ip, 0-1) < — (a - ip) > rns dp (24) 
da J -oo da 

where 0\ depends on the distance ~ \x e — x c \/9 (x e is the exit coordinate). Obviously the 
similar relations one can get for more complicated distribution of the initial angles. 
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Table 1: The initial conditions for calculations of the entering angle distributions 
and some results of these calculations 

R is the bending radius, l s is the total length of the scattering, < 6 > and < 56 > are the 
mean angle and standard deviation of the distribution. 

7 Calculations 

In this section we demonstrate applications of obtained relations for calculations of the 
specific cases volume reflection effect. All our calculations will be done for energy of the 
proton beam equal to 400 GeV, which was planned in CERN channeling experiment and for 
the (110) plane of the silicon single crystal. 

We start with detail consideration of calculation for proton beams, and then we also 
touch the case of negative particles. 

7.1 Atomic potential for calculations 

In contrast to channeling regime the above-barrier particle move throughout all the points 
of the potential and its influence on the motion of particle arises near critical point x c , where 
the particle is located close to top of potential barrier for long time. Because of this, we 
think that the correct choice of the one dimensional atomic potential is important problem. 

Taking into account aforesaid we choose the one dimensional atomic potential according 
to Ref.[16], where potential was calculated on the basis of the measurements of atomic form 
factors with the use of x-rays [17]. For this aim we use the expansion of the potential in the 
Fourier series. In a similar manner we find the one dimensional electric field and electron 
density. All these characteristics were calculated at room temperature. It should be noted 
that this potential was used for comparison of the theoretical description and experimental 
results for coherent bremsstrahlung of 10 GeV positrons for the (110) plane of the silicon 
single crystal [18] and in this work was shown very good agreement between theory and 
measurements. 

We also calculate the one dimensional potential and electric fields using the Moliere 
approximation of atomic form factors. The forms of the both potentials and electric fields 
are similar in between, but absolute values for the Moliere potential are larger approximately 
on 10 % than for experimental one. So, the potential bariers are equal to 23.36 eV and 21.38 
eV for Moliere and experimental approximations, correspondingly (see Fig. 2). 
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7.2 Calculation of volume reflection angle 

Fig. 3 illustrates the calculated volume reflection angle as a function of the transversal energy 
for several values of the radius. As expected, this function is periodic in accordance with 
Eq.(16). However, we see that a(E) is a discontinuous function. The purpose is obvious: it 
is a transition by bound from one potential cell to another. It is significant that function 
a(E) have a clear maximum and the value of the maximum decreases weakly with the radius 
decreasing, We find that the appearance of the maximum (or specifically, when derivative 
da/dE ^ 0) is connected with the form of the atomic potential. For some model potentials 
(such as a parabolic or polynomial [16] (with several numbers of terms) potentials) the 
maximum coincident with the point of discontinuity. For mentioned here model potentials is 
typical incorrect description of the electric field (electric field is a discontinues function near 
area of atomic location). Note that with the decreasing radius the interval of the possible 
angles increases. It means that the angle divergence is small at large enough radii. 

Knowledge of the function a(E) allow one to calculate the value a for given quantities 
x and 9. Fig. 4 illustrate the volume reflection angle as a function of the initial coordinate. 
Here the derivative da(x )/dx are also shown. There are two curves calculated at two 
different but very close in between initial angles 9 within angle interval which determined 
by Eq.(20). This fact shows the valuable variation of the a(a;o)-function at small variations 
of the initial angle. 

7.3 Initial distributions of particles 

Now we can find the initial distribution over the volume reflection angle and at the fixed 
initial angle with the help of Eq.(19). Fig. 5 illustrates our calculations. As follows from 
Fig. 4 there are three value of angle a where the derivative d9/dx is equal to zero (at fixed 
9). At these angles the distribution function dN/da tends to +oo. However, /°+ dN(ct) = 1, 
where a_ and a + are the minimal and maximal values of a. Comparing the two distribution 
at different values of 0-angle we see that only maximal value of a-angle is the same, but 
the couples of other angles (at which the derivative are equal to zero) are different. This 
fact is explanation of a disappearance of the two peaks at the next averaging over #-angle in 
accordance with Eq.(21). 

Our results of calculation according Eq.(19) we compared with the calculations by the 
Mote Carlo method of the same distributions and got very good agreement in between. 

The averaged over initial angles distributions are shown in Fig. 6 for several radii. As 
expected the peak at maximum angle is conserved. 

7.4 Entering angle distributions 

As was considered the multiple scattering distorts the initial particle distribution. For cal- 
culations of this factor we can find the distance where particle is scattered according to 
relation: l c = R9. Decreasing of angle 9 decreases the distance l c . On the other hand, the 
inequality l c ^> R9 C should be fulfilled. Thus, we can find that the optimal angle 9 lies in 
the angle range 510 -5 — 10~ 4 radian. It corresponds to l c = 0.125 — 0.25 cm for R = 25 
m. We also should take into account the distance l d after critical point (see Fig. 1). Our 
estimation (see Eqs.(23),(24)) show that the distribution form of the entering beam depends 
on the total distance l s — l c + Id and practically independent of the relation between l c and 
Id- Fig. 7 illustrates the results of calculations of the entering angle distributions at different 
radii. The initial (incoming) distribution we assume the physically narrow. The results of 
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calculations are partially reflected in the table 1. From our calculations one can see that the 
form of the outcoming angle distributions is close to Gaussian one. 

The angle width of the incident beam in CERN experiment is expected about some micro- 
radians. Then assuming its Gaussian form it is easy to find the resulting angle distribution 
after the bent single crystal. At the present time the precise Monte Carlo calculation of the 
angle distributions for conditions of the future experiment in CERN were made [19]. How- 
ever, we can say that our analytical description is consistent with these results and presents 
the additional information. 

7.5 Volume reflection of negatively charged particles 

We carried out some calculations of parameters of volume reflection effect for negatively 
charged particles at energy 400 GeV. In particular, we calculate the angle of volume reflection 
as a function of the transversal energy. For large enough radii ( pa 25m), the result is similar 
as in Fig. 3 but the value of maximal angle is less in pa 1.8 times. 

It should be noted, that one can expect some peculiarities in the volume reflection at 
small enough radii of bending for both sorts of charged particles. Really at very small radii 
the effect of volume reflection should be disappear due to dominant value of the linear term in 
effective potential. Mathematically, the next condition should be satisfied: SE ^> Uq, where 
Uq is the potential barrier of the usual one dimensional potential. For 400 GeV protons 
(antiprotons) the equality SE = Uq takes place at R ~ 3.5 m. We remind that through the 
paper only the (Oil) plane of silicon single crystal is considered. 

At the bending radii ~ 3.5 and less the angle of reflection (for both sorts of charged 
particles) may be of the different signs, or, in other words, the particles may be scattered 
in the opposite sides. Fig. 8 shows the angle of volume reflection as a function transversal 
energy for R = 1 m and negatively charged particles. One can see that the particles scattered 
mainly through relatively small angles, but there is the narrow range of a strong scattering 
through angles till 710~ 5 radian. Partially, we can explain this fact by the specific form of 
the effective potential (see Fig. 8b). We can also explain the discontinuity of a (E) -function 
(see Fig. 8a) by this reason. 

Now using Eqs.(23) and (24) we find the angle distribution after single crystal for physi- 
cally narrow initial beam. For calculations we select the thickness of the single crystal pa 0.01 
cm. Fig. 9 illustrates the results of our calculations. From here we get that pa 10, 3, 1 
percents of the beam scattered through angle more than 10~ 5 , 210~ 5 , 310 -5 radians, corre- 
spondingly. 

8 Conclusions 

Considered in paper analytical description of the process of volume reflection of charged 
particles in bent single crystals allow one to find necessary characteristic of scattered beam. 
We think that suggested method may be useful at different investigations at interactions 
particles with single crystals. This method may be used as additional to Monte Carlo 
calculations. Further development of our description is needed: 

1) the process of the volume capture of particles in channeling regime should be included 
in consideration with aim of the full description of the process; 

2) the extension of the method on single crystals with the variable radius is very desirable. 
We suppose that the first point may be easy enough realized because the main cause of 

the phenomenon is well-known. We believe that the second point may be also solved if take 
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into account that near the critical point the confined periodicity of the volume reflection 
angle takes place in some small volume of the bent single crystal. 



9 Acknowledgements 

Author would like to thank all the colleagues-participants of the CERN channeling experi- 
ment for attention and fruitful discussions. 

This work was partially supported by the grant INTAS-CERN 05-103-7525. 

References 

[1] A.M. Taratin, S.A. Vorobiev Physics Letter A, 425, (1987). 

[2] A.M. Taratin, S.A. Vorobiev Nucl. Instr. Meth. B47, 247 (1990). 

[3] R. P. Fliller et al, Phys. Rev. ST Accel. Beams 9, 013501 (2006). 

[4] N.V. Mokhov, et al., talk on channeling workshop, CERN, 8-9 March, 2006 

[5] Yu.M.Ivanov et al, preprint PNPI-2005-2649,Gatchina,2005. 

[6] W.Scandale et. al. First Observation of Proton Reflection from Bent Crystals. Talk on 
EPAC, Edinburgh, 26-30 June 2006. 

[7] CERN Channeling workshop, CERN, 8-9 March, 2006. 

[8] W. Scandale, In: Proceedings of the LHC Workshop. G. Gartskog and D. Rein Editors, 
Aachen, 1991, v.3, p.760. 

[9] A.D. Kovalenko, A.M. Taratin, E.N. Tsyganov . Preprint JINR-E1-92-9, Jan. 1992. 

[10] R.W. Assmann, S. Redaelli, W. Scandale. Optics Study for a Possible Crystal-based 
Collimation System for the LHC. Talk on EPAC, Edinburgh, 26-30 June 2006. 

[11] V. Guidi et. al. Characterization of Crystals for Steering of Protons through Channelling 
in Hadronic Accelerator. Talk on EPAC, Edinburgh, 26-30 June 2006. 

[12] M. Fiorini et. al. Experimental Study of Crystal Channeling at CERN-SPS for Beam- 
halo Cleaning. Talk on EPAC, Edinburgh, 26-30 June 2006. 

[13] Yu.A. Chesnokov, talk on channeling workshop, CERN, 8-9 March, 2006 

[14] E.N. Tsyganov. Fermilab Preprint TM-682, TM-684 Batavia, 1976. 

[15] X. Artru, S.P. Fomin, N.F. Shul'ga, K.A. Ispirian, N.K. Zhevago Physics Reports 412, 
89 (2005). 

[16] V.A. Maisheev, Nucl. Instr. Meth. B 119, 42, (1995). 

[17] D.T. Cromer and J.T. Weber, Acta Cryst. 18, 104, (1965); ibid. 19, 224, (1965). 

[18] N.K. Bulgakov et. al. Preprint JINR 1-83-621, Dubna,1983. 

[19] A. M. Taratin, talk on channeling workshop, CERN, 8-9 March, 2006 

9 



Figure 1: Scheme of the volume reflection of the proton beam. XYZ is the Carte- 
sian coordinate system at the entrance in single crystal, xyz is the local Cartesian 
coordinate system connected with the current location of the particle. V-axis is' 
directed normally to plane of figure. 9 and a are the initial and volume reflection 
angles. 
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Figure 2: One dimensional potential (a) and electric field (b) in the (110) plane 
of the silicon single crystal (at room temperature), as functions of the relative 
coordinate x/d, where d = 1.92 angstroms is the interplanar distance. The curves 
1 and 2 correspond to calculations based on the experimental and Moliere atomic 
form factors. 
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Figure 3: Angle of volume reflection as a function of the transversal energy for the 
values of radius equal to 25 m (a), 10 m (b) and 5 m (c), respectively. 
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a, rad 




Figure 4: Angle of volume reflection (a) and its derivative (b) as functions of the 
initial relative transversal coordinate x/d for radius equal to 25 m. The thick curves 
correspond to initial angle 6 = 1CT 4 radian and thin curves correspond to initial 
angle = 1CT 4 + 5.31CT 7 radian 
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Figure 5: Distribution of particles over the angle of volume reflection. Curves 1 and 
2 calculated for the two different initial angles (with the difference of these angles 
< 59) nearby 6 = 0.0001 radian and for R = 25 m. 
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Figure 6: Averaged over period of initial angles nearby 9 = 10~ 4 radian, distribu- 
tions of particles over the angle of volume reflection (a) . For visualization, the 
step distributions (b) were obtained by integration of the distributions (a) in angle 
intervals equal to 10~ 7 radian. Curves 1, 2 and 3 correspond to values of radius 
equal to 25, 15 and 10 m, respectively. For visualization, the y-coordinate for curves 
2 and 3 are increased in 5 and 10 times for (a) and in 5 and 12 times for (b) cases. 
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dN,1/rad 
da' 



a.rad 



Figure 7: Distributions of particles after volume reflection and multiple scattering 
(at the output of single crystal) as functions ( of the deflection angle. Zero angle 
corresponds to the initial direction of the beam. Calculation was made for different 
radii: 25 m (a), 10 m (b) and 5 m (c). For additional information see text and table. 
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Figure 8: One period of energy dependence of the volume reflection angle (a) and 
effective potential (b) for negative particles and for R = 1 m. The dotted line (in 
(a)) shows the discontinuity of a function. The horizontal lines (in (b)) correspond 
to energies at which the discontinuity takes place. 




a, radian 

Figure 9: Angle. distribution of the negative particles after bent single crystal (R= 
lm). The left part of the distribution are also shown with increasing of scale in 50 
times. The thickness of the single crystal is equal to 0.01 cm. 
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